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Abstract 

The spin-1 ferromagnetic Blume-Capel film with an alternating crystal field 
A = Ai on the odd layers and A = A2 on the even ones is considered in the mean 
field approximation. The ground state phase diagrams in the (di, c?2) plane (di = 
Ai/J and d2 = A2/J) are determined analytically; the number of their phases 
depends on the parity of the number of layers of the film. At finite temperature, 
fifteen types of topology, depending on the range of variation of di, are found 
in the {t,d2) plane for an even number of layers, but only fourteen for an odd 
number of layers. The phase diagrams exhibit a variety of multicritical points. In 
particular, a tricritical point C appears in the paramagnetic-ferromagnetic line 
of transition, but only for values of di larger than a threshold value d[^l. = 
which is but the tricritical crystal field (in the mean field approximation) of 
the spin-1 Blume-Capel model on a square lattice. Moreover, lines of transition 
presenting the reentrant and double reentrant behaviors may also take place. 
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1 Introduction 



The Blume-Capel model is a spin-1 Ising model with single-ion anisotropy. It has 
been originally introduced to study first-order magnetic phase transitions [1,2] 
and then applied to multicomponent fluids [3]. Later, it was generalized to the 
Blume-Emery-Griffiths model [4] to study the He^-He^ mixtures and a variety of 
other physical systems. 

The EC model has been investigated in detail using many approximate meth- 
ods, namely mean field approximation [1,2], high temperature series expansion 
[5], constant-coupling approximation [6], Monte-Carlo [7] and renormalization- 
group [8] techniques. All of these approximate schemes suggest the existence 
of a tricritical point at which the system changes from the second-order phase 
transition to the first-order one, when the transition temperature is plotted as a 
function of the parameter of anisotropy. 

On the other hand the EC model has been investigated on semi-infinite lat- 
tices with modified surface coupling. Several approximate methods have been 
used [9] (and references therein) . All of them show the possibility to have a phase 
with ordered surface and disordered bulk, and show the existence of the so-called 
'extraordinary', 'surface' and 'ordinary' transitions, and of multicritical points 
called special points [10]. More recently, with the development of the molecular 
beam epitaxy technique and its application to the growth of thin metallic films, 
renewed interest in thin film magnetism has been stimulated, and much exper- 
imental and theoretical efforts have been devoted to the subject [11]. The EC 
model with a rationally decreasing crystal field has recently been extended to a 
layered structure and analyzed in the mean field approximation [12]. The model 
exhibits a constant tricritical point and a reentrant phenomenon for a certain 
number of layers. 

Our aim in this work is to study, using mean field theory, the phase diagrams 
of a ferromagnetic spin-1 Blume-Capel film with an alternating crystal field: A = 
Ai on the odd layers and A = A2 on the even ones. The number of topologies 
of the phase diagrams is finite and depends only on the parity of the number of 
layers of the film. The phase diagrams present many multicritical points and may 
have transition lines exhibiting the reentrant and double-reentrant behaviors; this 
kind of behavior has been previously found in the Paramagnetic-Ferromagnetic 
transition line of the ordinary Blume-Capel model with a random crystal field 
[13] or with random nearest-neighbor interactions[14]. 
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The paper is organized as follows. In section 2 we present the model and 
give the mean field equations of the different order parameters. In section 3 we 
present the results, then we end in section 4 by a conclusion. 



2 Model and method 



The system we consider is the spin-1 Blume-Capel model of a film of N ferro- 
magnetically coupled layers with an alternating crystal field. The Hamiltonian 
of this system is given by 

n = -jY,SiS^ + Y.^i^i (1) 

(«j> » 

where Si = ±1,0; the first sum runs over all pairs of nearest-neighbors and Aj is 
the crystal field which takes the value Ai on the odd layers and A2 on the even 
ones. The mean field equations of state are straightforwardly obtained and are 
given by 

2sinh/?/ife 
= ^^'^ = e^A,+2cosh/3/., 

2\ 2cosh/?% 
- ^^'^ = e/^A.+2cosh/?^, 

hi = J{zmi + 1712) 

hk = J(mfe_i + zrrik + rrik+i), for 2 < A; < iV - 1 (4) 



with 



In these equations, (qk) is the reduced magnetization (quadrupolar mo- 
ment) of the layer; and z = 4 is the interlayer coordination number. The reduced 
total magnetization is m = J2k=i '^fc- The reduced free energy of the film is 
given by 




+-mi{Ami + 777.2) + -mN{mN-i + ^mN) (5) 
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3 Results and discussions 



The study of the phase diagrams in the (t,d2) plane exhibits the existence of 
different topologies depending on the range of variation of di {di = Ai/J, ^2 = 
A2/J and t = T/J are the values of the reduced crystal fields and temperature 
respectively). We note from the beginning that these topologies are limited in 
number and don't depend on the number of layers itself but only on its parity. 
In the following we will discuss in detail the case of an even number of layers 
and comment at the end on the odd number case. To understand the origin of 
the different topologies, we first determine analytically the ground state phase 
diagram in the (^1,^2) plane by looking for the lowest energy configurations. For 
an arbitrary even number of layers N = 2L, we distinguish six regions of this 
plane (Fig. 1(a)) (but only five for an odd number of layers (Fig. 1(b))): 

(a) di < 3, d2 < S and di + ^2 < 6 — 1/L; the stable configuration^ is (11)^, 
where the spins of the sites of all the layers are all equal to 1. 

(b) 3 < di < 4, ^2 < di and (L — l)di + Ld2 < 6L — 4; the stable configuration 
is 01(11)^"^, where the spins of the first layer are all equal to while those 
of the other layers are all equal to 1. 

(c) di > 4 and ^2 < 2; the stable configuration is (01)^ where the spins of the 
odd layers arc all equal to while those of the even ones are all equal to 1. 

(d) 3 < d2 < 4:, di < d2 and Ldi + (L — 1)^2 < 6L — 4; the stable configuration 
is (11)^-^10 . 

(e) di <2 and d2 > 4; the stable configuration is (10)-^. 

(f) di > 2, d2 > 2, di + d2 > 6 - 1/L, [L - l)di + Ld2 > 6L - 4 and 
Ldi + {L — l)d2 > 6L — 4; the stable configuration is the non magnetic state 
(00)^. 

From Fig. 1(a) it is clear that for fixed values of di, the system exhibits various 
transitions by varying d2 ■ For the following ranges of variation of di : i) < cZi < 
2, ii) 2 < di < 3 - 1/L, iii) 3 - 1/L < di < 3, iv) 3 < di < 4 and v) di > 4, the 
system presents the following phase transitions respectively: 

i) from the state (11)^ to the state (11)^^^ 10 at ^2 = 3 and then to the state 
(00)^ at d2 = 4. 

^Wc use obvious notations for the model where the natural basic entity of the different states and 
phases is a bilayer 
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ii) from the state (11)^ to the state (11)^ "^10 at (^2 = 3 and then to the state 
(00)^ at Ld2 = (6L - 4) - (L - l)di. 

iii) from the state (11)^ to the state (00)^ at ^2 = (6 - l/L) - di. 

iv) from the state 01(11)^"^ to the state (00)^ at (L - 1)^2 = (6L - 4) - Ldi. 

v) from the state (01)^ to the state (00)^ at ^2 = 2. 

It is worth to note here that, contrary to the two or three dimensional Blume- 
Capel cases (or the case of a film with a homogenous crystal field), the ground 
state of certain layers may be magnetic even for relatively large values of the 
crystal field in these layers. For example, for di < 2, the even layers (except the 
final one) are in the magnetic states m2k = 1 even with a crystal field, in these 
layers, reaching the value d2 = 4: (which is larger than the critical crystal field 

(2) ('Vt 

for a square {dc =2) or three dimensional cubic (dc =3) lattices). 

For finite temperatures the transition lines are determined by solving numer- 
ically the equations [J] which may have more than one solution; the one which 
minimizes the free energy |21 corresponds to the stable phase. Transitions of the 
second order are characterized by a continuous vanishing of the magnetizations 
while those of the first order exhibit discontinuities at the transition points. 

In order to describe the different entities in the phase diagrams, the Griffiths 
notations [15] will be adopted: the critical end-point B'^A^ denotes the inter- 
section of m lines of second order and n of first order; the multicritical point C 
denotes the intersection of m lines of second order; the tricritical point, which 
is the intersection of a line of second order and a line of first order, is denoted 
in particular by C . In addition, the following notations of the different phases 
will be adopted: (FF)^ is, for example, a state where all the layers are in the 
ferromagnetic state, i.e. the magnetization and the quadrupolar moment of each 
layer are both different from zero. {FF)^~^FO is a phase where all the layers are 
in the ferromagnetic state except the last one which is in the non magnetic state 
(noted by 0) and is characterized by vanishing magnetization and quadrupolar 
moment. (PP)^ is a phase where all the layers are in the paramagnetic state i.e. 
the magnetization of each layer is equal to zero but not the quadrupolar moment 

For an even number of layers {N = 2L), fifteen types of topology of the phase 
diagram in the {t,d2) plane are found. They may be gathered into six groups, 
each one is characterized by a certain range of variation of di . 
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1) For < di < (e.g. Fig.2(a)), there is a critical line separating the 
paramagnetic phase (PP)^ from the ferromagnetic ones {FF)^, {FF)^~^FO 
and (FO)^. These ferromagnetic phases are separated by two transition lines of 
first-order at low temperatures which meet two critical lines at two tricritical 
points (totally inside the ordered phases) Ci and C2 respectively. The first- 
order lines reach the ^2 axis at the fixed values (independent of L) 1^2 = 3 and 
(^2 = 4 respectively. The critical lines present both a reentrant behavior near 
the paramagnetic-ferromagnetic (P-F) transition line and meet this last at two 
multicritical points (-B^)i and (.6^)2 respectively. The topology described here 
is characterized, contrary to the ordinary Blume-Capel model, by the absence of 
a tricritical point in the (P-F) transition line . The value of di over which this 
point appears is found to be d[fl « 1.848 « which is but the mean field 
tricritical crystal field of the spin-1 Blume-Capel model on a square lattice. 

(2) 

2) For df/i^ < di < 2, four types of topology are found with increasing di, 
all of them present a tricritical point in the P-F transition line and present 
the same features of the type-1 topology with some differences. For di in the 
vicinity of df^]. the lines of transition (F-F lines) separating the ferromagnetic 
phases (FF)^, {FF)^''^FO and (FO)^ still exhibit the reentrant phenomenon 
near the P-F transition line and meet this one at the multicritical points (5^)i 
and (-8^)2 (e.g. Fig2(b)). With increasing di, the F-F lines meet the P-F line 
at the multicritical point and at the critical end-point BA^ respectively, and 
still present the reentrant behavior (e.g. Fig. 2(c)); for higher values of di, the 
former lines lose the reentrant behavior and meet the P-F line at two distinct 
critical end points {BA'^)i and {BA^)2 as is depicted in Fig. 2(d). Finally, for di 
close to 2 the P-F first-order line develops a reentrant topology (e.g. Fig.2(e)). 

Note that for di < 2 the ferromagnetic phase (FO)^ is stable for an infinitely 

large value of d2 at temperatures lower than a temperature tc (di) which is 

(2) 

independent of the number of layers 2L. tc {di) is but the critical temperature 
of the ordinary two-dimensional spin-1 Blume-Capel model corresponding to di. 
This may easily be understood even beyond the mean field approximation. In 
fact at low temperatures and for large values of d2, the odd layers are all non 
magnetic {m2k+i = 0) and the system of equations of the magnetizations m2k 
of the even layers is decoupled to L equivalent two-dimensional equations with 
a crystal field di; the magnetizations m2k (and then the total magnetization) 
vanish at the same temperature which is tc (di). For di > 2 the ferromagnetic 
phase is unstable for large values of d2, as will be seen, in accordance with the 
ground-state phase diagram. 
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3) For 2 < di <3 — 1/L, three types of topology exist. For these topologies 
the ferromagnetic state (FO)^ disappears from the ordered phases. The ferro- 
magnetic phases {FF)^ and {FF)^~^FO are, as above, separated by a transition 
line of first-order, at low temperatures, which reaches the d2 axis at ^2 = 3 
and meets the second order line at a tricritical point Ci. This last line meets 
the P-F line at a critical end-point BA^. For di in the close vicinity of 2, the 
P-F line of transition exhibits a reentrant and double-reentrant behaviors (e.g. 
Fig.2(f)); with increasing di, only the double reentrant still exist (e.g. Fig.2(g)) 
and then it disappears for higher values of di (e.g. Fig.2(h)). Note that the 
P-F transition line reaches, for 2 < di < 3 — 1/L, the d2 axis at d2 given by 
(L — l)d2 = (6L — 4) — Ldi; the surface of the {FF)^~^FO phase then decreases 
with increasing di and it collapses at di = 3 — 1/L. 

4) For 3 — 1/L < di < 3, the phase diagram is of the same type as that of the 
ordinary spin-1 Blume-Capel model (e.g. Fig.2(i)); the first-order line reaches 
the ^2 axis at c?2 given by ^2 = (6 — 1/L) — di . 

5) For 3 < di < 4, three types of topology are found. These are characterized 
by the appearance of the ferromagnetic phase OF{FF)^'~^ at low temperatures. 
The ferromagnetic phases (FF)^ and OF[FF)^^^ are separated by a critical 
line of approximately constant temperature which meets the P-F transition line 
at a critical end-point BA^. The P-F line reaches the d2 axis at d2 given by 
Ld2 = (6L — 4) — (L — l)di . For di far from 4, a typical phase diagram is given in 
Fig.2(j). With increasing di, the P-F transition line develops a double-reentrant 
behavior (e.g. Fig.2(k)), and when di approaches 4 this line presents a reentrant 
and double-reentrant behaviors (e.g. Fig2(l)). 

6) For di > 4, three types of topology are fond; these are characterized by 
the appearance, in addition to the state OF{FF)^-'^, of the state (OF)^ at low 
temperatures. The ferromagnetic phases (FF)^, OF{FF)^~^ and (OF)''" are sep- 
arated by two critical lines of approximately constant temperatures which meet 
the P-F transition line at two critical end points {BA'^)i and {BA'^)2 respec- 
tively. The P-F line reaches the d2 axis at d2 = 2. For values of di very close 
to 4 the P-F line presents a reentrancc and a doTibIc reentrancc (e.g. Fig.2(m)); 
the reentrancc disappears with increasing di (e.g. Fig.2(n)) and then the double 
reentrancc disappears for high values of di (e.g. Fig.2(o)). 

For an odd number, iV = 2L -|- 1, of layers, the ground-state (Fig.l(b)) has 
only five different phases (instead of six for N even); in fact, for small values of di 
{di < 2) only two phases now exist: (ll)-^l (for d2 < 4) and (10)^1 (for d2 > 4). 
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For finite temperatures, fourteen types of topology exist in this case: 

i) For < di < dj^^ the topology is like the type-1 above (in the case of an 
even number of layers) but with only two ferromagnetic phases {FF)^F 
and {FO)^F separated by a unique line having a tricritical point Ci at low 
temperatures. 

ii) For d[^]. < di < 2 there are four types of topology which are like those of 
the group 2 above, but with the modification mentioned in i). 

iii) For 2 < di < 3, there arc no topologies like those of the group 3 above. The 
topologies, in the number of three, are all of the ordinary Blume-Capel type, 
but for the two first ones {di close to 2) the P-F transition line develops a 
reentrance for di very close to 2, and a reentrance and double-reentrance 
for higher values of di (Fig. 3). 

iv) For 3 < di < 4, there are three types of topology which are like those of 
the group 5 above, but with the ferromagnetic phase OF[FF)^~^0 at low 
temperatures (instead of OF(FF)^~^ for N even). 

v) For di > 4, there are three topologies which are like those of the group 6 
above, but with the ferromagnetic phases (OF)^O and OF{FF)^~^0 at low 
temperatures (instead of (OF)^ and OF{FF)^~^ for N even). 



4 Conclusion 



We have studied, using mean field theory, the phase diagrams of a ferromagnetic 
spin-1 Blume-Capel film with an alternating crystal field: A = Ai on the odd 
layers and A = A2 on the even ones. The ground state phase diagrams in the 
{di,d2) plane (di = Ai/J and d2 = A2/J) arc determined analytically; the 
number of their phases depends on the parity of the number of layers of the 
film. At finite temperatures, fifteen types of topology, depending on the range 
of variation of di , are found in the {t,d2) plane for an even number of layers, 
but only fourteen for an odd number of layers. The phase diagrams exhibit 
many multicritical points. In particular, a tricritical point C appears in the 
paramagnetic-ferromagnetic line of transition, but only for values of di larger 
than a threshold value d[^]. = which is but the mean field tricritical crystal 
field of the spin-1 Blume-Capel model on a square lattice. Moreover, lines of 
transition presenting the reentrant and double-reentrant behaviors also appear 
in the phase diagrams. 
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Figure captions 



Fig.l The ground state phase diagrams in the {di,d2) plane for an arbitrary 
number N of layers, (a) N = 2L, {h) N = 2L + 1. 

Fig. 2 The phase diagrams in the (t,c?2) plane for a generic even value of N 
{N = 6) and generic values of di. (a) di = 1., (b)di = 1.85, (c) di = 1.855, 
(d) di = 1.9, (e) di = 1.99, (f) di = 2.01, (g) di = 2.1, (h) di = 2.25, 
(i) di = 2.75, (j) di = 3.5, (k) di = 3.8, (1) di = 3.9, (m) di = 4.01, (n) 
di =4.1, (o) di = 5. The solid lines are of second-order, the doted ones are 
of first-order. 

Fig.3 The phase diagram in the (t,d2) plane for a generic odd value of N (iV = 7) 
and generic values of di in the range 2 < di < 3. The number accompanying 
each curve is the value of di. The solid lines are of second-order, the doted 
ones are of first-order. 
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Fig.2 
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